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The chiral Fe-based langasites represent model systems of triangle-based frustrated magnets with 
a strong potential for multiferroicity. We report neutron scattering measurements for the multichiral 
Ba3MFe3Si20i4 {M = Nb, Ta) langasites revealing new important features of the magnetic order 
of these systems: the bunching of the helical modulation along the c-axis and the in-plane distortion 
of the 120° Fe-spin arrangement. We discuss these subtle features in terms of the microscopic spin 
Hamiltonian, and provide the link to the magnetically-induced electric polarization observed in these 
systems. Thus, our findings put the multiferroicity of this attractive family of materials on solid 
ground. 

PACS numbers: 75.85.+t, 75.25.-j, 75.30.Ds, 75.30.Gw 


The quest for novel magnetic orders that promote mul¬ 
tiferroicity and strong magnetoelectric coupling drives a 
considerable amount of current research (see e.g. [1-4] 
for recent reviews). In this respect, Fe-based langasites 
are propitious materials displaying an original “multi¬ 
chiral” magnetic order in which a 120° Fe-spin arrange¬ 
ment undergoes an additional helical modulation [5-8]. 
In fact, the magnetoelectric effect in these systems has 
been demonstrated at both static [9] and dynamical levels 
[10]. In addition, Zhou et al [11] have reported the emer¬ 
gence of ferroelectricity at the onset of magnetic order. 
This ferroelectricity has subsequently been confirmed by 
Lee et al. [12] , who have also demonstrated the magnetic- 
field manipulation of the electric polarization. However, 
these works report two different ferroelectric axis at zero- 
field and, above all, leave the precise mechanism promot¬ 
ing such a ferroelectricity unidentified. 

Magnetically-induced ferroelectricity in many multifer- 
roic cycloidal magnets is due to the spin-orbit coupling 
(SOC) via the so-called inverse Dzyaloshinskii-Moryia 
(DM) mechanism [13-15]. However, this mechanism is 
ineffective for the reported helical magnetic structure 
of the langasites [16, 17]. An electric polarization ob¬ 
served in magnetic helices rather than cycloids has been 
theorized [18], as well as observed experimentally in 
RbFe(Mo 04)2 for instance [16, 17, 19]. However, the in¬ 
voked mechanisms still rely on the SOC and appear inap¬ 
plicable in the langasite case. Bulaevskii et al. [3, 20, 21] 
have discussed another fundamental mechanism that can 
generically operate in triangle-based frustrated magnets. 
The proposed mechanism is independent of SOC as it 


simply relies, in the language of Mott insulators, on high- 
order hopping-term contributions that enable the redis¬ 
tribution of the electronic charge. This purely electronic 
mechanism can be supplemented by a contribution from 
lattice degrees of freedom [3]. In any case, it requires a 
distortion of the perfect 120°-spin structure and, despite 
its generality, has not been verified experimentally so far. 

In this paper, we carefully analyze neutron scatter¬ 
ing experiments to revise the magnetic structure of the 
Ba 3 MFe 3 Si 2 0i4 (M = Nb, Ta) langasites. We discover 
additional features indicating the bunching of the he- 
licoidal modulation and significant deviations from the 
120°-spin structure in these systems. In order to explain 
these features, we propose an extended spin Hamiltonian 
incorporating different magnetic anisotropies. These ex¬ 
tra terms prove to be crucial for understanding the 
magnetically-induced ferroelectricity in the Fe langasites. 
Furthermore, we find that this ferroelectricity emerges 
from a magnetoelectric coupling that is equivalent to the 
Bulaevskii et al. proposal at the phenomenological level. 

The Fe langasites crystallize in the non-polar non- 
centrosymmetric space group P321, in which the mag¬ 
netic Fe^+ ions {S = 5/2) form a triangular lattice of tri¬ 
angles in the a6-plane (see Fig. 1). The magnetic order 
of these systems, observed below the Neel temperature 
Tat, has been rationalized within the spin Hamiltonian 
[6, 7, 22-24]: 
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Here the first term takes into account the exchange 
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also measured on IN22 at 2 K with a final wavevector 
fixed to 2.66 and at zero energy transfer in order to 
get a quantitative determination of the elastic structure 
factors of the magnetic satellites. The crystal was ori¬ 
ented with the a zone axis vertical so as to probe the ( 6 *, 
c*) scattering plane. Two types of measurements were 
carried out using different monochromator-analyser con¬ 
figurations: Heusler-Heusler (H-H) and Graphite-Heusler 
(G-H). The former allows us to discriminate the magnetic 
and chiral contributions from the nuclear one, whereas 
the latter increases the incident neutron flux and allows 
us to determine the chiral scattering from an unpolarized 
beam [7], hence with a larger fiux. 


FIG. 1: (Color online) (a) Ba 3 NbFe 3 Si 20 i 4 crystallographic 
structure in the a5-plane and along the c-axis, with five mag¬ 
netic exchange paths Ji to J 5 . (b) Schematic representation 
of the revisited magnetic structure projected along the c-axis. 
The 3-fold axis is broken leading to different Ji interactions 
(case of an isosceles triangle represented). The bunching of 
the helices in the Nb compound is reported through the suc¬ 
cessive angles between consecutive magnetic moments along 
c, deviating from the regular angle 27rT=51.4°. 


interactions, where the sum runs over the spin pairs con¬ 
nected by different paths shown in Fig. 1 (a). The 
nearest- and next-nearest-neighbor interactions Ji and 
J 2 in the a 6 -plane are both antiferromagnetic. This gen¬ 
erates magnetic frustration within the Fe triangles and 
causes the 120 ° spin structure. Among the three anti¬ 
ferromagnetic exchange couplings between adjacent ab- 
planes, the strongest one is J 3 or J 5 depending on the 
structural chirality. These couplings generate the heli¬ 
cal modulation of the 120 ° structure along the c-axis. 
The propagation vector of this modulation depends on 
the balance between the inter-plane interactions J 3 , J 4 
and J 5 . The second term in Eq. (1), where the sum 
runs over intra-triangle spins, corresponds to the addi¬ 
tional DM interactions with the DM vector (D) along 
the c-axis, which is allowed by symmetry. This DM term 
restricts the spins to the a 6 -plane and favors one intra¬ 
triangle chirality as observed experimentally [22, 23, 25]. 
Alternatively, a weak easy-axis single-ion anisotropy can 
have similar effect [26]. 

We have revised this structure by carrying out single¬ 
crystal neutron scattering measurements on two mem¬ 
bers of the Fe langasite family with Nb or Ta on the 
la Wyckoff site. The two single-crystals, grown by 
floating-zone method in an image furnace [9] , display op¬ 
posite structural chiralities [25]. The experiments were 
done on two spectrometers installed at the Institut Lane 
Langevin: IN5 time-of-flight spectrometer using unpolar¬ 
ized neutrons, and GRG-GEA IN 22 triple-axis spectrom¬ 
eter, equipped with CRYOPAD, using polarized neutrons 
and allowing longitudinal polarization analysis. The dy¬ 
namical structure factors were measured on IN5 at 1.5 K, 
for both crystals in rotation about the a zone axis, with 
an incident wavelength of 4 A. The Nb compound was 


The intensity maps at zero energy cut measured on IN5 
for the Nb and Ta compounds are shown in Eigs. 2(a) 
and (b) respectively. The Neel temperature of these com¬ 
pounds are nearly identical (27 K for Nb and 28 K for Ta), 
while the propagation vectors are r = (0,0,0.1429) and 
(0,0,0.1385) respectively [9, 27]. The latter is obtained 
from the position of the magnetic satellites at Q = \H^r\ 
in Eigs. 2 . In addition to these first-order satellites, we 
clearly observe higher-order satellites at Q =| ^ ± 2t | 
and Q =\ ^ ± 3r I in both compounds suggesting that 
the helical modulation is distorted. 

Polarized neutron measurements with longitudinal po¬ 
larization analysis on the Nb compound allowed us to 
extract the nuclear (ctat), the magnetic (cfm) and the 
chiral {(Jch) scattering [25], hence to identify the nuclear 
and/or magnetic nature of these extra peaks. This anal¬ 
ysis reveals that the second-order satellites are mainly of 
nuclear origin. This is usually attributed to magnetoelas¬ 
tic effects and will not be further discussed. In contrast, 
the magnetic contribution is dominant in the third-order 
satellites [see Eig. 2(c)], which confirms the deformation 
of the magnetic helices. 

This deformation allows further explanation of the sub¬ 
tle details in the spin-wave spectrum. As we can see in 
Eig. 3, the spectra of the Nb and Ta compounds display 
two branches emerging from the magnetic satellites [ 22 ] . 
The high energy branch corresponds to spin components 
in the a 6 -plane while the low-energy branch is associated 
to spin components along the c-axis. The main features 
of this spectrum were initially captured by the model 
of Eq. 1 using the standard Holstein-Primakov formal¬ 
ism in the linear approximation (see Eig. 4 of ref. 22 ). 
However, this simple model fails to reproduce a spectral 
weight extinction observed in the a 6 -branch around 2.3 
meV in both compounds, and a tiny extinction around 
1.7 meV in the c-branch, visible only in the Ta compound 
(red and orange arrows respectively in Eig. 3) [25]. 

In order to reproduce the deformation of the he¬ 
lices leading to all these features, we have considered a 
straightforward extension of Eq. 1 . We included extra 
single-ion anisotropy terms with respect to a local 2 -fold 
zioc axes, i.e. along a, b and —a — b for the three ions of 
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FIG. 2 : (color online) Zero energy cut (integrated from -0.1 to 0.1 meV) from IN5 measurements at 1.5 K, on the Nb (a) and Ta 
(b) compounds. The blue arrows on the energy maps indicated the [0 ±1 ^]* directions of the Q-ID cuts shown below in vertical 
log scale, (c)-(d) summarize the measurements performed on the Nb compound at 2 K on IN22 with polarized neutrons, (c) 
shows the magnetic (red), nuclear (green) and chiral (blue) elastic contributions from scans obtained by rotating the crystal 
(ct;-scans) with the H-H configuration. The solid lines are gaussian fits. The magnetic and chiral contributions measured in 
the H-H (circles) and G-H (triangles) configurations for the main (0, 1, i ± r) and extra (0, 1, i ± 3t) satellites, versus the 
calculated contributions using the model given in the caption of Fig. 3, are displayed in (d). Note that the G-H configuration 
cannot separate the magnetic from nuclear scattering, but we checked that ctm ctn ^ cfm for these satellites. 



[ 0 - 1 ^] [0 1 ^] 

FIG. 3: (Golor online) Spin-waves spectra of Nb (a) and 
Ta (b) measured along the [0 -1 ^]* and [0 1 ^]* directions 
respectively at 1.5 K on IN5. Below: Spin-waves calcu¬ 
lations with the following parameters for the Nb(c)/Ta(d): 
Ji = 0.85/0.85 ± 0.1 meV, J 2 = 0.24/0.24 ± 0.05 meV, 
J 3 = 0.053/0.22 ± 0.03 meV, J 4 = 0.017/0.014 ± 0.05 meV 
and J 5 = 0.24/0.053 ± 0.05 meV, -0.028/-0.028 ±0.005 
meV, easy-plane anisotropy K — 0.045/0.053 ± 0.005 meV. 
The red/orange arrows show the spectral weight extinctions. 


a triangle respectively: 

= ( 2 ) 


We first determine the magnetic ground state of the ex¬ 
tended Hamiltonian by means of a real-space mean-field 
energy minimization. Then we compute the correspond¬ 
ing magnetic and chiral elastic scatterings for the first 
and third order magnetic satellites, and finally we com¬ 
pute the corresponding spin-wave spectrum [25]. 

A quantitive determination of the parameters of the 
revised Hamiltonian can be achieved for the Nb com¬ 
pound. The intensity of the the first-order (0, k, ^ ± r) 
and extra (0, k, ± 3r) satellites (Figs. 2 (d)) increases 
with the single-ion anisotropy and agrees with the ex¬ 
periments (see Fig. 2) for |iF|=0.045±0.05 meV. This 
term is also at the origin of the anti-crossings observed 
in the spin-waves (at 2.3 meV for the a6-branches and 
at 1.7 meV for the c-branches): these spectral weight 
extinctions occur at the crossing point of branches with 
the same symmetry, emerging from the ±r and — r elas¬ 
tic magnetic satellites, which are mixed by the single ion 
anisotropy KS‘1 [25, 28]. The anti-crossings are more 
pronounced in the Ta case (see Fig. 3 (b)) which qual¬ 
itatively agrees with a higher value of K [25]. Another 
signature of the presence of tiny anisotropies is a gap in 
the c-branch, which reaches its energy minimum at ~ 
0.4 meV for Nb [26] and ^ 0.2 meV for Ta as shown 
in Fig. 4 (a). The variation of this gap is actually op¬ 
posite for the easy-plane {K > 0) and for the easy-axis 
{K < 0) anisotropies. In the latter case, the c-gap in¬ 
creases with increasing \K\ and is already around 1 meV 
for iF=0.045 meV, which is larger than the values ob¬ 
served experimentally. These values are actually com¬ 
patible with an easy-plane anisotropy combined with a 
DM term (with D < 0): the c-gap increases when \D\ 
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increases and decreases when K increases. The Nb c-gap 
is reproduced for T)=-0.028 meV and if=0.045 meV (see 
Fig. 4(b)). The smaller gap of the Ta compound is in 
turn compatible with a larger K than in the Nb one (see 
Fig. 4(c)). This analysis yields the extended model given 
in the caption of Fig. 3 with a set of parameters repro¬ 
ducing the spin-waves for the Nb and Ta compounds (see 
Figs. 3(c-d)). Moreover, in the Nb case these parameters 
are fully compatible with the magnetic structure factors 
and with the Electron Spin Resonance measurements [23] 
that have been reanalyzed carefully [25]. In real space, 
this revised model globally corresponds to a deviation 
of the angle between adjacent spins along the helix axis 
away from the ideal 27rr = 51.4° value: It gets smaller 
near the local easy-planes and larger far from them (see 
Fig. 1 (b)). Thus, we conclude that the magnetic struc¬ 
ture of the Fe langasites corresponds to a bunched helix, 
as reported in some rare-earth intermetallics [29]. 

Next, we analyze in more detail the first-order satel¬ 
lites along the [0 0 ^]* line. These (0, 0, ^ ± r) satel¬ 
lites are clearly visible in our data (see Fig. 2), where 
they display the same temperature dependence as the 
(0, 1, ^ ± r) magnetic ones [25]. In fact, the longitudi¬ 
nal polarization analysis confirms their magnetic origin 
(see Fig. 2(c)). The observation of these satellites is very 
surprising since, in the case of a perfect 120 ° spin ar¬ 
rangement within the Fe triangles, there should be an 
extinction of the magnetic structure factor along the [0 
0 ^]* line. We note, that the forbidden (0, 0, r) satel¬ 
lite has also been observed by magnetic resonant X-ray 
diffraction (together with 2r and 3r satellites in the first 
Brillouin zone) [30] . This was attributed to a modulated 
out-of plane (butterfly) component of the magnetization 
produced by the in-plane intra-triangle DM vector. How¬ 
ever, neutrons are blind to this component along [0 0 
^]* since they are only sensitive to the magnetic compo¬ 
nents perpendicular to the scattering vector Q. Thus, 
the observation of these ( 0 , 0 , satellites in our neu¬ 

tron experiments clearly demonstrates that the relative 
orientation of the magnetic moments within the trian¬ 
gles deviates from 120 °, hence suggesting that the intra¬ 
triangle Ji interactions are not equivalent. Our neutron 
data cannot single out whether the distortion of this tri¬ 
angular arrangement is local (helically modulated from 
planes to planes) or global. However, only the latter 
case, implying a loss of the 3-fold symmetry axis of the 
P321 structure, is compatible with a macroscopic polar¬ 
ization measured along the two-fold axis. The loss of the 
3-fold axis was already invoked from Mossbauer experi¬ 
ments and x-ray magnetic diffraction [27, 30], and was an 
important ingredient for explaining the magnetoelectric 
excitations revealed by means of THz spectroscopy [10]. 
In addition, it is compatible with the phonon spectrum 
[31]. This finding has strong implications on the multi- 
ferroic properties of these systems, as we discuss in the 
following. 

The P321 space group of the Fe langasites allows the 
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FIG. 4: (Color online) (a) Neutron scattering intensity versus 
the energy at the magnetic satellite position for Ta and Nb, 
obtained from ID Q-cut (integration ranges A/i = A/c = A/ = 
0.05) measured on IN5 with a wavelength of 8 A. Below: Cal¬ 
culated scattering intensity of the c-branch only showing the 
c-gap of the Nb (b) and Ta (c) compounds using the revised 
model with the parameters of the caption of Fig. 3. 


magnetoelectric coupling of the form: 


Fme = a{P,{S^-S2yS^- 



[(Si+S2)-S3-2Si.S2]}. 


where y and x stand for the a-axis and its perpendicu¬ 
lar direction in the aPplane, and Si to S 3 are the three 
spins on a triangle (see Fig. 1). This coupling can have 
a purely electronic origin as discussed in [ 20 ] or emerge 
from magnetostriction effects in the symmetric Heisen¬ 
berg exchange [3]. In any case, at the phenomenolog¬ 
ical level, its form is identical to the one discussed in 
[3, 20, 21]. The presence of such a coupling means that 
a spin-induced electric polarization has to be expected 
whenever the spin arrangement deviates from the perfect 
120° structure [32]. This polarization will lie in the ab- 
plane. Thus, the new magnetic structure deduced from 
our experiments is compatible with the polarisation mea¬ 
surements of Lee et al [ 12 ] rather than those of Zhou et 
al. [11]. This is also in tune with the DFT-based calcu¬ 
lations reported in [33] , in which a magnetically-induced 
polarization was also obtained in the a 6 -plane due to the 
reduction of symmetry of the magnetic structure to C 2 ’. 
In this work, however, this reduction was the result of 
spin-orbit coupling (likely via symmetric anisotropic ex¬ 
change) , and the precise magnetoelectric coupling behind 
such an electric polarization was not identified. 

The enhancement of the electric polarization in the 
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a6-plane by the application of a magnetic field has also 
been reported in [12]. This can be understood on the 
basis of the above coupling, although its connection with 
the revised magnetic structure is more subtle. In this 
case, the 120° spin-structure of each triangle is addition¬ 
ally deformed by the Zeeman coupling, which induces a 
(weak) ferromagnetic component and thus an extra local 
contribution to the electric polarization. In the perfect 
helical structure, however, these individual field-induced 
electric polarizations coming from the different layers of 
the helix rotate along the c axis in such a way that they 
cancel each other. In the bunched helical structure, in 
contrast, this rotation is modified so that the cancela¬ 
tion is avoided, which eventually enables an overall finite 
polarization induced by the magnetic field. This will hap¬ 
pen even if the initial (zero-field) structure corresponds 
to a perfect 120° arrangement within the Fe triangles. 
Thus, the magnetic-field induced electric polarization es¬ 
sentially relies on the bunching of the helix along the 
c-axis. 

By carrying out neutron scattering experiments, we 
have identified subtle features in the magnetic order of Fe 
langasites, with drastic consequences on their properties. 
The first feature is a bunching of the helical modulation 


related to single-ion magnetocrystalline anisotropy, and 
the second one is the distortion of the 120° spin structure. 
At the phenomenological level, we have established that 
the magnetically-induced electric polarization at zero and 
finite magnetic fields results from a generic magnetoelec¬ 
tric coupling that is equivalent to that put forward by 
Bulaevskii et al Our study thus sets the basis for under¬ 
standing the multiferroic properties of Fe langasites rep¬ 
resenting model-case magnetically-induced ferroelectrics 
in triangle-based frustrated magnets. Finally, the close 
relation between the multifaceted chiral properties of the 
Fe langasites and their multiferroic potential is expected 
to motivate further studies. 
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VALIDATION OF MAGNETIC ANISOTROPY TERMS BY ELECTRON SPIN RESONANCE 

Since our present neutron study has resulted in substantially different magnetic-anisotropy parameters of the inves¬ 
tigated Ee langasites compared to our previous study that was based on the electron spin resonance (ESR) [1], it is 
essential to check the compatibility between these two complementary techniques. The currently derived parameters 
are roughly an order of magnitude larger than the previous estimates. It should be stressed though that the previous 
ESR analysis only took into account both sources of magnetic anisotropy separately. As our current study has shown, 
the easy-plane single-ion anisotropy {K > 0) effectively reduces the zero-field gap of the c-branch that is opened 
by the out-of-plane c-component of the DM interaction, D, which thus explains the reduced anisotropy parameters 
proposed in Ref. 1. 



B(T) 


FIG. 1: The frequency-field diagram of measured (symbols) and simulated (lines) AFMR modes in Ba3NbFe3So20i4 at 4 K. 
The simulations correspond to the c-component DM anisotropy constants D = —0.028 meV and the easy-plane single-ion 
anisotropy K = 0.045 meV. The gap of the c-branch determined by neutron measurements is indicated. 









2 



T(K) 


FIG. 2: Temperature dependence of the magnetic contributions of the (0, -1, l-r) and (0, 0, r) satellites measured by neutron 
scattering and rescaled for comparison. 


We have repeated the calculation of the AFMR modes detected in the previous ESR experiment [1]. As can be 
seen in Fig.l, with no adjustable parameters an excellent agreement is found for the values D = —0.028 meV and 
K = 0.045 meV, which convincingly validates the results of our neutron scattering study. 

We stress that the in-plane a6-component of the DM interaction, Dab, does not affect the lowest AFMR modes and, 
therefore, can not be determined from the AFMR experiment. However, it does affect the ESR line-width anisotropy 
in the paramagnetic phase, which was found to saturate at the value of = 1.27 at high temperatures [1]. 

The above minimal model of the magnetic anisotropy, employing only D and K yields AH^^/AHc = 0.75, therefore 
an additional anisotropy term is needed to explain the experiment. The isotropic a term of the single-ion anisotropy 
cannot boost this ratio above 1, while \Dab\ = 2.5|D| reproduces the experiment. The minimal model of the magnetic 
anisotropy in Ba3NbEe3Si20i4 that reproduces all the experiments performed so far thus includes D = —0.028 meV, 
\Dab\ = 0.07 meV and K = 0.045 meV. 


DEVIATION FROM A 120° MAGNETIC ARRANGEMENT AND OTHER MODELS 

The magnetic structure of the Ee langasite is described by three magnetic helices propagating perpendicular to the 
(a, b) plane with a propagation vector (0, 0, r), based on magnetic moments at 120° from each other in the a6-plane. 
The magnetic structure factor of this structure, at the reciprocal space position H = ±r with H a reciprocal lattice 
vector, writes : 


Fm{Q = H±t)=p T] fi\Q\)'Fi±f,ueiQ.ft, 

u=l,3 


with p=0.2696 10 cm, / the magnetic form factor of the Ee^+ atoms, fv the position of the atom u in the cell: 
fi = (x,0, ^), r 2 = (0,x, and fs = (—x, —x, with x=0.2496. The Eourier components of the magnetization 
distribution for the atom u is with fi the magnetic moment amplitude, {u^v) and orthonormal 

basis in the a6-plane, and (j)j, a phase which should account for the relative 120° dephasing between the three Bravais 
lattices. 

Eor H = (0,0,^), the magnetic structure factor becomes zero: 


Fm(0,0,^±t) =pf{\Q\) 


E 

u=l,3 


fiu ± ifiv 


0^'^^V 




= pf{\Q\)[ 


l2U±il2V 


(1 






Eor in-plane magnetic moments, the only way to observed a magnetic signal on the (0, 0, ^±r) satellites, implies a 
deviation from the 120° arrangement of the magnetic moments. The measured temperature dependence of the (0, 0, 
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r) extra-satellites, produced by the deformation of the 120° magnetic arrangement is shown in Fig. 2. It follows the 
temperature dependence of strong first order magnetic satellites. 

For completeness, we have tried two other models: the first one involves, in addition to the Hamiltonian of Eq.l, a 
biquadratic term of the form [{Si.S 2 ){S 2 — Ss).Si + S 2 -Ss{Ss — Si).S 2 + {Ss.Si){Si — 5'2)-5'3] due to the spin-phonon 
coupling, (ii) the second model, called KSEA (for Kaplan-Shekhtman-Entin-Wohlman-Aharony), takes into account 
terms quadratic in which leads to experimental observations such as 3r satellites and spectral weight extinction in 
the spin-waves [4, 5]. We found that none of these models quantitatively reproduce our data, therefore we discarded 
them. 


NOTE ON THE CHIRALITY 

The excellent agreement between the calculated and observed spectral weight extinction in the spin-waves when 
including a single-ion anisotropy in the Hamiltonian is only observed for one triangular chirality (anticlockwise 120 ° 
rotation of the spins on each triangle), in agreement with previous determination of the magnetic chirality [3]. Both 
Ta and Nb single-crystals thus have the same triangular chirality. Their structural chirality is opposite, as determined 
from anomalous X-ray scattering. This structural chirality, in turn, imposes the helical chirality to be opposite in 
both compounds. 


NEUTRON POLARIMETRY 


The simplest neutron polarimetry technique to probe the chiral scattering was first introduced by Moon, Riste 
and Koehler [2], and is called longitudinal polarization analysis (LPA). In this case the final and initial polarizations 
are parallel, which can be achieved typically on a triple-axis spectrometer with polarizing monochromator/analyser 
{e.g. Heusler crystals). In addition, two flippers select the polarization states + and — (parallel or antiparallel to the 
polarization axis). The spin-flip terms (scattering processes changing the sign of the polarization). 
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and the non-spin-flip terms (scattering processes leaving the sign of the polarization unchanged). 
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of the partial differential cross-section can be measured independently. 

In neutron polarimetry, a right-handed coordinated system is usually chosen with the x—axis along the scattering 
vector Q, the axis in the scattering plane and the 2 :—axis perpendicular to the scattering plane, so that the magnetic 
interaction vector has zero x component. In the following, we shall use the simplified notations: the nuclear cross- 
section cttv is proportional to \N\‘^ with N the nuclear structure factor. The magnetic cross-section ctm is proportional 
to with M±_ the magnetic interaction vector, which is the projection of the magnetic structure factor Em onto 

the plane perpendicular to the scattering vector Q. The chiral cross-section dch is proportional to Mch defined as: 


Mch = i{Ml^Ml - 


These contributions can then be easily determined from linear combinations of the spin-flip and non-spin-flip cross- 
sections, knowing the initial polarisation Pi : 
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Mch oc 


2 R 


Alternatively, one can use an initial unpolarized beam (produced by a graphite monochromator for instance) and 
perform polarization analysis. In this case, the chiral scattering can be derived from the difference: 


. 0 - 


t 0 + 


= Mch 


where the symbol 0 is meant for zero beam polarization of the incoming neutrons. 

The flipping ratio of R ~ 16 was determined experimentally by measuring the neutron counts in the Spin-Flip (SF) 
and Non Spin-Flip (NSF) channels on different nuclear Bragg peaks. This value was taken into account in the IN22 
data analysis to correct the intensities measured in the SF and NSF channels from the imperfectly polarized neutron 
beam. 


DETAILS ON THE NEUTRON DATA ANALYSIS 

The analysis of the neutron data was performed using the spinwave software developped at LLB by S. Petit. The 
starting point is a bilinear spin Hamiltonian 

= '^CEF + 2 ^ 

where J'ij denotes the effective exchange coupling (including standard exchange, anisotropic, Dzyaloshinskiy-Moryia 
or bond directed coupling); ^ single ion anisotropy term : 

Hc^P = 3B2o {Sl-S{S + l)). 

At the mean field level : 

'Hmf = 

i 

3 

Starting from a random configuration for the observables (-Si), the contribution to 1-Lmf at site i (in the unit cell) 
is diagonalized in the Hilbert space defined by the {\Sz)} ,Sz = —S',..., S basis vectors (2S + 1 vectors). This yields 
the energies Ei^n and the wave functions |0i,n)- The updated expectation values, (S^)', at each step of the iteration 
procedure, are given by: 




o — Ei^nlksT 


■ (0i,n |Si |0i 


with 


= E ksT) . 

n 

These are then used to proceed to site j, and this is repeated until convergence. 

In the present case, the starting point is a spin arrangement with a magnetic cell 7 times larger than the nuclear 
one, which is a correct description for the Nb langasite and only approximate for the Ta langasite (propagation vectors 
respectively equal to ~l/7 and ^1/7.22). This allowed us to find that, without single-ion anisotropy, the mean field 
procedure converges to a perfect helix where the moments are at 120° from each other. Including the single ion 
anisotropy we find that the helix becomes bunched (still with a 120° magnetic order within the triangles). Then, 
imposing the loss of the 3-fold axis (by taking different Ji on one triangle), we find that the 120° arrangement is 
perturbed. After this mean field step, the same software computes the neutron magnetic structure factors for all 
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the satellites that can be compared with the measured ones. Last the spin-waves spectrum (based on the Holstein- 
Primakov description) of the minimized spin configuration is obtained and compared to the experiment. We found 
an excellent agreement between the spin configuration imposed by our model Hamiltonian, the extra satellites and 
the observed extinction in the spin-waves. The single-ion anisotropy produces a bunching of the helix, the rise of 3r 
satellites, and the extinction in the spin-waves. The loss of the 3-fold axis additionally produces a deformation of the 
120° arrangement and the rise of the first order satellites along the [0 0^]* direction. 

The biquadratic coupling 


V — ^ijkl 

ijkl 


Si • Sj Sk • Si 


can be taken into account in the software. Assuming the symmetry properties: 


^ijkl ^jikl 
^ijkl — ^ijlk 
^ijkl ~ ^klij 

the non-linear nature of the biquadratic coupling is treated by a decoupling scheme, ending up with bilinear couplings 
only. The following approximation is performed: 

( V ^ Y^ij ^ab ^kl ^ijkl Vk Vl ^ab + 4 ^ikjl Vk Vl)) Sj 

\m = {Si) 

leading, from the formal point of view, to a new effective anisotropic exchange term. 


[1] A. Zorko, M. Pregelj, A. Potocnik, J. van Tol, A. Ozarowski, V. Simonet, P. Lejay, S. Petit, and R. Ballou, Phys. Rev. Lett. 
107, 257203 (2011). 

[2] R. M. Moon, T. Riste and W. C. Koehler, Phys. Rev. B 181, (1969) 920. 

[3] V. Simonet, M. Loire, and R. Ballou, Eur. Phys. J. Special Topics 213, 5 (2012). 

[4] L. Shekhtman, Amnon Aharony, and O. Entin-Wohlman et al, Phys. Rev. B 47, (1993) 174. 

[5] A. Zheludev, S. Maslov, G. Shirane et al, Phys. Rev. B 59, (1999) 11432. 



